An analysis of /o — cr mixing in light cone QCD sum rules 



A. Gokalp *, Y. S and O. Yilmaz ' 

Physics Department, Middle East Technical University, 06531 Ankara, Turkey 

(February 2, 2008) 



Abstract 



We investigate fo — cr mixing in the framework of light cone QCD sum rules 
and by employing the experimental results about the decay widths and the 
masses of these scalar mesons we estimate the scalar mixing angle by using a 
meson data and fo meson data. The two values we thus obtain of the scalar 
mixing angle are not entirely consistent with each other, possibly indicating 
that the structure of these mesons cannot be simple quark-antiquark states. 
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The hght scalar mesons have been the subject of continuous interest in hadron spec- 
troscopy. In addition to well established isoscalar /o(980) and isovcctor ao(980) scalar 
mesons, the recent experiments and theoretical analysis find evidence for the existence of an 
isoscalar (t(600) [1], and an isodoublet ft;(900) scalar mesons [2]. Although the structure of 
these light scalar mesons have not been unambiguously determined yet [3], the possibility 
may be suggested that these nine scalar mesons below 1 GeV may form a scalar SU(3) flavor 
nonet [4] . In this case, one would expect a mixing in the scalar sector similar to that occurs 
in the pseudoscalar sector. One thus expresses the physical (t(600) and /o(980) states as 
linear combinations of isoscalar ss and {uu + dd) / \/2 states as 

^ UU + dd . ^ 

a — cos Os 7= sm Og ss 

v2 

uu + dd _ 
fo = sm6'5 — --j= h cos6's ss , (1) 

where Og is the scalar mixing angle. 

In this work, we estimate the scalar mixing angle Og by employing light cone QCD sum 
rules. The method has been used for the calculation of hadronic coupling constants and for 
the studies of hadronic properties [5] . We study the hadronic decays fo — > tttt and a ^ nn 
within the framework of this method and by utilizing the experimental data about these 
decays and about the masses and the widths of /o and a mesons we deduce the scalar mixing 
angle Og. 

In order to study 5'7r7r-vertex, where S denotes the scalar meson fo or a, we consider the 
two point correlation function 

T,{p,q) ^ijd'x e^^-^ < 7r+{q)\T{j;-{x)f{0)}\0 > , (2) 

where p and are the four-momentum and the interpolating current for 7T~ meson, is 
the interpolating current for scalar fo or a meson, and q is the four-momentum of tt"*" state. 
The interpolating currents in terms of quark fields are the axial vector (x) — It(x)7^75d(x) 
for pion, and the scalar currents j^''{x) = (s'mOg/ \/2)[u{x)u{x) + d{x)d,{x)] + cos^^ s{x)s{x) 
for fo meson and j"'{x) = (cos Og / ^/2)[u(x)u(x) + d{x)d(x)] — sin^^g s(x)s(x) for a meson. 
The scalar current for fo and a meson is assumed to have a no n- vanishing matrix element 
between the vacuum and a scalar meson state, < S \ \ >= Xs, where Xs is called the 
overlap amplitude. The overlap amplitude A^. for a meson [6] and the overlap amplitude A/^ 
for fo meson [7] were previously determined by using two-point QCD sum rules method. 

The correlation function T^(p, g) can be written in terms of two independent invariant 
functions Ti and T2 by using Lorentz invariance as 

T^(p, q) = iT, ((p + p^) p^ + {{p + q)\ p^) q^ . (3) 

We consider the invariant function Ti((p + q)^, p^). In order to construct the theoretical 
part of the sum rule in terms of QCD degrees of freedom we calculate the function Ti by 
evaluating the correlation function in the deep Euclidean region, where [p -\- qY and p^ 
are large and negative, as an expansion near the light cone = 0. This expansion involves 
matrix elements of non-local quark-gluon operators between pion and vacuum states which 
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defines pion distribution ampUtudes of increasing twist. We retain terms up to twist four 
accuracy since higher twist amplitudes are suppressed by powers of (p+g)^] or (p)^] 
and thus arc known to make a small contribution [8]. In our calculation we use the full light 
quark propagator with both perturbative and nonperturbative contributions which is given 
as [9] 

iS{x, 0) = < 0\T{q{x)q{0)}\0 > 
. 4^ <qq> 



X 



^2n^x^ 12 



Y^-o <m> 

dul4ai,,G'"'{ux) - 4:iu^G^''{ux)jA + ... . 



x^ 



(4) 



We note that the two-point correlation function T^, and therefore the invariant function Ti, 
have contributions coming only from the first terms of the interpolating currents for /o or 
(7 mesons. In other words, since the pion does not have any strangeness content only the 
{sm9s/V2){uu + dd) part of the interpolating current for /q meson and the {cos 9s /V^){uu + 
dd) part of the interpolating current for cr meson make any contribution in the evaluation 
of time-ordered products of pion and scalar meson currents in Eq. 2. We, therefore, in our 
calculation of the theoretical part of the sum rule consider the scalar current as j^{x) = 
\u{x)u{x) + d{x)d{x)]/y/2, and denote the resulting expression for the invariant function Ti 
as T[. After a straightforward computation and after performing the Fourier transforms the 
invariant function T[ is obtained to twist four accuracy as 



+ 
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Ml 



1 I {p + uq)-q 



3 (p + uq)^ 
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\p + q{l-ai- va3)Y 



Ml 



\p + q{ai + vas)] 




(5) 



In the above expression the functions (fianiu) and (ppniu) are the twist 3 quark- antiquark 
pion distribution amplitudes which are defined by the matrix elements [10] 



< 



7:+{q)\u{x)i-f5d{0)\0>= due''"i-''^pAu) , 



(6) 



< n^{q)\u{x)a,,^,d{m >= (l - ^) (Q,^^ ' ^m?-) £ «"''>,,(k) , (7) 

where /i^r = M^/2mq with rriq = {ruu + Tnd)/2 is the twist 3 distribution amplitude nor- 
malization factor. Although in the evaluation of the light quark propagator given in Eq. 
4 we put iriu = rrid = 0, we like to note that niq in the normalization factor is ob- 
tained by making use of the Gell-Mann-Oakes-Renner relation [11] as nig — —f^M^/4{qq). 
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We work in the fixed-point gauge x'^A^ — 0, consequently the path-ordered gauge factors 
P exp ^igs Jq da Xfj,A'^(ax)^ which appear in between the quark fields and which assure 
gauge invariance are not included in the matrix elements [9]. The twist 3 quark-antiquark- 
gluon pion distribution amplitude (p3Tr{ai) is defined as [10] 

X j L>o;i(/P3.(aOe'^"^"^+""^^ (8) 

where Dai = daida2da3S{l — ai — a2 — a^) . After performing the Borel transformation with 
respect to the variables Ql = —{p + qY and Ql = —p^, we obtain the theoretical expression 
for the invariant function T[ in the form 
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(9) 



where and M| are the Borel parameters and 
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The correlation function Tj^(p, g") satisfies a dispersion relation, therefore we can represent 
the invariant function Ti as 



T,{{p + q)\ / jdsds' 



p^""\s, s') 



{p + q)'^)]{s' - p^ 



(10) 



where the hadronic spectral density p'*"'^(s, s') receives contributions from the single-particle 
states and also from higher resonances and continuum states. We therefore saturate this 
dispersion relation by inserting a complete set of one hadron-states into the correlation 
function and we consider the single-particle tt and scalar meson S, where S denotes /o or cr, 
states. This way we obtain the expression for the invariant function Ti in the form 



T^{{p + q)\ p') 



< I J'm I ^^(^) I ^ ^< S{p + q)\j'^ I > 



+ ds ds' 



so 



[{p + q)^-Ml] {p^-MD 
p'^\s,s') 



[s-{p + q)^]{s'-p^) ' 



(11) 



where now the hadronic spectral density p^'"^{s, s') includes the contributions of higher 
resonances and the hadronic continuum which contribute in a domain T> of the {s, s') plane 
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starting from two thresholds sq and s'q. The matrix element < tt'^tt \ S > defines the 
coupling constant gsT^n 

< 7r+(g)7r-(p) | S{p + q) >= gs^n (12) 
and the current-particle matrix elements are given as 

<Sip + q)\f \ 0>^Xs , (13) 



<0\j;' \^^{p)>=^Up, ■ (14) 

The contribution coming from the continuum can be identified by using global quark-hadron 
duality [12] with the QCD contribution above the thresholds Sq and s[y This way the pole 
contribution in which the coupling constant g^TrTr appears is isolated. After performing the 
same Borel transformation that was applied to the theoretical expression for the invariant 
function Tl we obtain for the hadronic representation of the same invariant function the 
result 

T.iMl M|) = Xs U gs.. e"^!/^? e"^^/^! 

+ I ds I rfs'p'=°"*(s,s') e-^/*^? e-^'/^2 . (15) 

■'SO J s'q 

The sum rule for the coupling constant g57r7r then follows in accordance with the QCD sum 
rules method strategy by equating the expressions Ti(M^,M|) obtained for the invariant 
function Ti((p + qY^p^) by theoretical calculation given in Eq. 9 and by physical consider- 
ations which is given in Eq. 15. For this purpose we note that the expression for T[ given 
in Eq. 9 should be multiplied by cos Og to obtain the sum rule for g^^^^ and by sin Og for the 
corresponding sum rule for g/oTTTr- In order to do this we have to identify the second term 
in Eq. 15 representing the contribution of the continuum with a part of the term calculated 
theoretically in QCD. We follow the prescription, for this so called the subtraction of the 
continuum, for the cases where the Borel parameters corresponding to channels with differ- 
ent mass scales cannot be constrained to be equal [13,14]. This prescription is based on the 
observation that the distribution amplitudes '{^^^{u) and (Paixiu) are polynomials in (1 — m), 
therefore in order to compute their contribution in the duality region T> we can write 

vvM + 1^ = J: bk{i - uf , <p,M - ^ ^ = E - uf . (16) 

k — k — 

Since the contribution of the twist 3 quark-antiquark-gluon term in Eq. 9 is small, we thus 
affect the continuum subtraction in the leading twist 3 quark-antiquark term. Therefore, we 
finally obtain the sum rules for the coupling constants go-TTTr and g/o7r7r in the form 

ga^n ^cosOg g'^^^ , gf^^^ ^ sin g'f^^^ (17) 

where 
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(18) 



where A — Sq/M"^ with Sq the smallest continuum threshold. 

In the numerical evaluation of the sum rule we use the twist 3 pion distribution ampli- 
tudes given by [10] 
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(20) 



</'37r(Q;i) = 360aia2a3 



l + ^(7as-3) 



(21) 



where C!^{2u — 1) are the Gcgcnbaucr polynomials. The overlap amplitudes Xs were deter- 
mined previously by employing two-point QCD sum rules method as = (0.17±0.02) GeV^ 
[6] and Xf^ — (0.18 ± 0.02) GeV'^ [7]. We also adopt the values of the parameters in the 
distribution amplitudes at the rcnormalization scale 1 GeV as f3Tt{lGeV) = 0.0035 GeV^, 
iOsTri^GeV) = —2.88, and 03 = [10]. Moreover, we also use the values {qq){lGeV) = 
-(0.240 GeVf and U = 0.132 GeV [13]. The mass of the a meson is taken as 
Mg. = (483 ±31) MeV [1] in which statistical and systematic errors are added in quadrature. 

We then study the dependencies of the sum rules for the couphng constants g'^^^ and 
g'fgT^T^ on the Borcl parameters Mf and M| by considering independent variations of these 
parameters. As for the continuum threshold parameter sq we use the value sq = 1.2 GeV^ 
for a meson [6] and Sq = 1.1 GeV^ for Jq meson [7] that were obtained in the calculation 
of the overlap amplitudes A^. and A/^ using two-point QCD sum rules method. We find 
that the sum rules are quite stable for the variation of the Borel parameters in the ranges 
2 < Mf < 6 GeV'^ and 0.4 < M| < 0.8 GeV'^. The variation of the couphng constant 
g^^^ as a function of Borel parameters Mf and M| is shown in Fig. 1 and that of g/gT^Tr 
shown in Fig. 2. The values of the coupling constants thus can be obtained by choosing 
the value of Mf about at the middle of the stability region as 3.2 < g'^^^ < 3.9 GeV and 
3.4 < g%,^ < 3.8 GeV. 

On the other hand, the decay widths of a and /o mesons can be calculated in terms of 
the coupling constants ga-T^Tv and g/QTTTr which are defined through the matrix elements given 
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in Eq. 12. The decay width S — > tt+tt can thus be obtained in terms of the couphng 

constant g57r7r ^ 



where p* = y M| — (M^+ + Mt^-Y /2 is the magnitude of the three-momentum of the either 
of the pions in the CM. frame. The experimental width 

V^'iS ^ tttt) = ^ T{S ^ TT+TT-) = 3 r(5 ^ 7r°7r°) (23) 

is given as T/,, = 40 - 100 McV with M/„ = (980 ± 10) MeV [15] for /o meson, and a best 
fit to the Dahtz plot of — > Ti'^a — > tt+tt+tt" decay in E791 Collaboration experiment [1] 
results in = (483 ±31) MeV and T„ = (338 ± 48) MeV for a meson. If we use these 
experimental values we then obtain the coupling constants g57r7r as ga-T^Tr — (2-6 ± 0.2) GeV 
and g/oTTTT = (1.6 ±0.8) GeV. The scalar mixing angle 6s can then be obtained by using Eq. 
17 and utilizing on the one hand the values of the coupling constants gsinr determined from 
experimental results and on the other hand the values of g'^^^ estimated by light cone QCD 
sum rules method. This way we determine the scalar mixing angle Og from a meson results 
as 9s — (41 ± 11)° and from /o meson results as 9s — (27 ± 13)°. These two values do not 
exclude each other, since they have a small overlap region, but given the large difference in 
the average values of the scalar mixing angle 6s determined using cr meson and /o meson 
properties we may assert that our results indicate that the structure of a and /o mesons 
cannot be simple quark-antiquark states. 

It should be mentioned that the previous determinations of the scalar mixing angle 9s 
through experimental data coming from charmed meson decays also lead to a somewhat 
inconsistent picture. The J/ip /o(98O)0 and J/ip /o(980)cj decays were used to 
estimate the mixing angle [16] and the result 9s = (34 ±6)° was obtained, which is consistent 
with our results. The analysis of the experimental results Dg — > /o(980)7r"'" and Dg — > 07r+ 
[17] on the other hand yields the mixing angle in the range 35° < 9s < 55°. In these analysis, 
it is the mixing angle in /o meson state that is determined. However, such a mixing angle 
implies that a meson also has a strange component, therefore the decay Dg a7r~^ should 
also occur, but this decay is not observed to make any contribution to the Dg —>■ tt+tt+tt" 
decay by the E791 Collaboration [18]. Therefore, this apparent inconsistency about the 
scalar mixing angle may be taken to indicate that the structure of scalar mesons a and /o 
are more complicated than simple quark-antiquark states, a conclusion that is also supported 
by our calculation. 
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FIG. 1. The coupling constant g^^^ as a function of the Borel parameter for different values 
of the Borel parameter M|. The curves denote the limits of the stability region. 
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FIG. 2. The coupling constant gj^TTTr ^® ^ function of the Borel parameter for different 
values of the Borel parameter Ml- The curves denote the limits of the stability region. 
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